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ABSTRACT 


As  part  of  a  general  investigation  of  two- stream  super¬ 
sonic  diffusers,  a  theory  was  developed  for  the  bounded  tur¬ 
bulent  mixing  between  coaxial  high  velocity  streams  of 
different  composition  and  temperature.  The  usual  assumption 
of  negligible  axial  pressure  gradient  is  invalid  in  this  case 
because  the  outer  stream  is  bounded  by  a  cylindrical  duct.  A 
highly  simplified  flow  model  is  used  which  retains  the  essential 
features  of  the  actual  mixing  process;  the  resulting  system  of 
equations  is  amenable  to  manual  calculation. 

The  theory  is  correlated  with  data  obtained  from  an 
annular  nozzle  wind  tunnel  having  a  central  supersonic  core  of 
rocket  exhaust  gases.  Correlation  is  also  made  with  Mikhail's 
low  velocity  data  for  which  hie  more  elaborate  theory  is  avail¬ 
able  for  comparison.  The  simplified  theory  predicts  with 
reasonable  accuracy  the  experimental  axial  static  pressure  and 
velocity  distributions. 


3 


AEDC-TR-61-18 


COKTENTS 

Page 

ABSTRACT .  5 

NOMENCLATURE .  7 

INTRODUCTION .  9 

DEVELOPMENT  OF  THEORY .  9 

EXPERIMENTAL  INVESTIGATION .  17 

DISCUSSION  OF  RESULTS .  18 

CONCLUSIONS .  21 

REFERENCES .  21 

APPENDIX  -  Iteration  Procedure  for  Mixing  Analysis  ....  23 

TABLE 

Tunnel  Dimensions  and  Operating  Parameters .  25 

ILLUSTRATIONS 

Figure 

1.  Schematic  of  Tunnel  for  Two-Stream  Diffuser 

Investigation .  27 

2.  Schematic  of  Velocity  Profiles  for  Coaxial 

Stream  Duet  Flow  , ' .  28 

3.  Approximation  of  Actual  Velocity  Profiles  for 

Coaxial  Duct  Flow .  29 

4.  Mixing  Model  for  Calculation  of  Local  Rates 

of  Heat  and  Momentum  Transfer .  30 

5.  Dividing  Streamline  Velocity  Ratio .  31 

6.  Mach  Number  Function  for  Air  and  Typical  Rocket 

Products .  32 

7.  Variation  of  Mixing  Zone  Spreading  Rate  with 

Mach  Number  for  Zero  Secondary  Velocity .  33 

8.  Variation  of  Turbulent  Mixing  Parameter  with  Mach 
Number  for  Fully  Developed  Mixing  of  an  Axlsym- 

metric  Jet .  34 


5 


AEDC-TR-61-18 


Figure  Page 


9.  Schematic  of  Mikhail's  Experimental  Apparatus  ....  35 

10.  Comparison  of  Theory  with  Mikhail's  Experi¬ 
mental  Data .  36 

11.  Comparison  of  Theoretical  and  Experimental 

Centerline  Velocity  Distribution .  37 

12.  Schematic  of  Supersonic  Jet  Expanding  into 

Supersonic  Stream .  . .  38 

13.  Axial  Static  Pressure  Distribution  for  Mach 
Number  3.  07  Tunnel  at  Toa  =  1100“R 

a.  Poa  40.  8  psia . 39 

b.  Poa  ~  23,  5  psia .  40 

c.  Poa  °  16.  0  psia .  41 

14.  Theoretical  Pressure  and  Mach  Number  Distri¬ 
butions  for  Mach  Number  3.  07  Tunnel . .  42 

15.  Theoretical  Outer  Stream  Mach  Number  Distri¬ 
bution  for  Mach  Number  3.  07  Tunnel .  43 


6 


AEDC-TR-61-1B 


A 

b 

Cfw 

'P 

D 

F 

J 

k 

M 

P 

Po 

€ 

R 

To 

u 

w 

X 

Y 

y 

€ 

P 

a 

T 


NOMENCLATURE 

Cross -sectional  area 

Width  of  mixing  zone 

Skin  friction  coefficient  at  wall 

Specific  heat 

Diameter 

Total  stream  thrust 

Joule's  constant,  mechanical  equivalent  of  heat 

Turbulent  mixing  parameter 

Mach  number 

Static  pressure 

Total  pressure 

Total  energy  transfer  rate 

Local  energy  flux  between  streams 

Gas  constant 

Total  temperature 

Axial  velocity 

Weight  flow 

Axial  distance  from  initiation  of  mixing 
Coordinate  in  mixing  model 
Ratio  of  specific  heats 
Apparent  kinematic  viscosity 
Density 

Turbulent  mixing  parameter 
Shear  stress 


SUBSCRIPTS 

1 .  2  Refers  to  cross  section 

a  Outer  stream  conditions 
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Arithmetic  average  conditione  over  an  interval 

Conditions  along  dividing  streamline 

Duct 

Conditions  in  core  flow 

Maximum  velocity  adjacent  to  mixing  region 
Minimum  velocity  adjacent  to  mixing  region 
Conditions  at  duct  wall 
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INTRODUCTION 


An  experimental  and  theoretical  investigation  of  supersonic  two- 
stream  diffuser  operation  is  currently  in  progress  at  the  Rocket  Test 
Facility  {RTF),  Arnold  Engineering  Development  Center  (AEDC),  Air 
Force  Systems  Command  (AFSC).  The  investigation  is  concerned  with 
operation  of  an  axisymmetric  wind  tunnel  having  a  central  core  of  rocket 
exhaust  gases.  Such  a  diffuser  is  normally  operated  with  a  fully  devel¬ 
oped  duct  shock  system  which  causes  rapid  mixing  of  the  streams.  Over¬ 
all  diffuser  performance  may  be  accurately  predicted  for  such  a  case 
by  using  one- dimensional  theory. 

At  certain  operating  conditions,  however,  it  is  desirable  to  main¬ 
tain  the  back-pressure- induced  shock  system  far  downstream  in  the 
diffuser,  and  for  this  case  the  mixing  characteristics  of  the  two  streams 
while  both  streams  remain  supersonic  become  important.  Simple  one¬ 
dimensional  theory  is  inadequate  for  analysis  of  the  supersonic  mixing 
process,  and  since  rigorous  analysis  of  the  bounded  two-stream  mixing 
problem  is  very  difficult,  a  simplified  theory  for  bounded  mixing  of 
coaxial  compressible  streams  has  been  developed  and  is  presented  in 
this  report.  A  later  report  will  be  presented  on  the  general  two- stream 
diffuser  investigation.  The  theoretical  and  experimental  results  on 
two- stream  mixing  are  also  applicable  to  other  coaxial  stream  mixing 
problems  -  for  example,  problems  concerning  the  operation  of  jet  pumps 
and  propulsion  system  thrust  augmentation  devices. 

The  simplified  theory  was  correlated  with  experimental  data  obtained 
at  AEDC  using  the  tunnel  configuration  shown  in  Fig.  1  operating  at  low 
back  pressure  and  with  low  velocity  two- stream  mixing  data  obtained  from 
other  agencies. 


DEVELOPMENT  OF  THEORY 


Generally,  one  of  two  approaches  is  taken  to  solve  free  turbulent 
mixing  problems.  The  first  approach  is  to  write  equations  of  motion  for 
the  flow  field  and  find  solutions  satisfying  the  bovindary  conditions.  The 
second  approach  is  to  assume  velocity  profile  shapes  and  then  use  integral 
methods  to  define  the  flow  field.  No  purely  theoretical  explanation  for 
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turbulent  processes  has  yet  been  devised,  and  at  least  one  empirical 
constant,  a  turbulent  mixing  parameter,  must  be  used  for  either 
approach. 

Por  the  problem  of  interest,  the  mixing  of  two  streams  inside  a 
duct,  the  usual  assumption  of  constant  pressure  mixing  is  invalid. 

Mikhail  (Ref.  1),  in  his  treatment  of  incompressible  mixing  of  coaxial 
streams  in  a  duct,  is  one  of  the  few  investigators  to  consider  axial 
pressure  gradients.  In  his  case,  in  spite  of  the  simplifying  condition  of 
incompressible  flow,  a  number  of  assumptions  were  necessary  to  solve 
the  integral  equations,  and  three  empirical  turbulent  coefficients  were 
used  to  correlate  theory  with  experimental  data.  In  addition  to  axial 
pressure  gradients,  it  is  necessary  to  include  effects  of  compressibility, 
energy  transfer  between  the  streams,  wall  friction,  and  different  stream 
composition  in  a  realistic  treatment  of  the  behavior  of  a  rocket  jet  ex¬ 
hausting  into  a  bounded  coaxial  supersonic  airstream.  Considering 
Mikhail's  difficulties  with  a  much  less  complex  case,  it  is  not  surprising 
that  the  inclusion  of  these  effects  in  the  more  classical  approaches  to 
turbulent  mixing  problems  results  in  equations  which,  if  they  can  be 
written,  are  extremely  difficult  to  solve,  even  with  numerical  methods. 

In  view  of  these  difficulties,  an  approach  was  taken  during  the  present 
study,  which  is  a  considerable  departure  from  earlier  attempts  to  solve 
such  nonconstant  pressure  mixing  problems. 

A  schematic  of  the  actual  flow  model  for  coaxial  stream  duct  flow 
is  shown  in  Fig.  2.  The  axial  variation  of  flow  conditions  in  this  duct 
may  be  treated  approximately  by  assuming  that  the  mass  of  each  stream 
remains  separated  and  that  each  is  one- dimensional  at  any  section  (Fig.  3). 
If  relations  for  the  shear  stress  and  heat  transfer  between  the  streams 
and  at  the  duct  wall  may  be  written  in  terms  of  the  one- dimensional 
stream  properties,  then  one- dimensional  equations  for  momentum,  con¬ 
tinuity,  and  energy  may  be  written  for  each  stream  using  the  unknowns, 
Toa»  Toj,  p,  Aj,  Ma,  and  Mj. 

The  assumptions  of  the  analysis  are: 

1.  Flow  is  steady. 

2.  Local  rates  of  heat  and  momentum  transfer  between  the  core  and 
outer  flow  are  the  same  as  for  mixing  between  respectively 
uniform  two-dimensional  streams. 

3.  Each  stream  is  initially  one- dimensional  and  remains  one- 
dimensional  with  no  subsequent  mass  interchange. 

4.  Each  stream  acts  as  a  perfect  gas  having  constant  specific  heat 
and  composition. 
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5.  Turbulent  shear  stress  over  a  finite  interval,  Ax,  is  the  average 
of  the  shear  stress  at  the  end  points  of  the  interval. 

6.  The  turbulent  Prandtl  number  for  the  mixing  process  is  unity. 

7.  The  velocity,  temperature,  and  composition  profiles  across  the 
mixing  zone  are  linear  (Pig.  4). 

8.  The  dividing  streamline  separating  the  mass  flow  of  the  respec¬ 
tive  streams  occurs  at  the  midpoint  of  the  mixing  zone  velocity 
profile. 

9.  The  turbulent  exchange  parameter,  k,  is  constant  for  the  entire 
mixing  process. 

10.  The  static  pressure  is  constant  across  any  duct  cross  section. 

1 1 .  The  flow  at  the  duct  wall  is  turbulent,  and  the  wall  friction  coef¬ 
ficient,  cf„ ,  is  constant  over  the  entire  duct. 


TURBULENT  TRANSFER  QUANTITIES 


The  local  turbulent  exchange  in  the  mixing  zone  was  calculated 
using  the  simple  mixing  model  of  Fig.  4  for  mixing  between  uniform 
streams.  The  shear  stress  in  the  turbulent  region  is 


where  ^  is  the  apparent  turbulent  kinematic  viscosity.  For  turbulent 
Prandtl  number  unity,  the  apparent  viscosity  for  energy  transfer  is  iden¬ 
tical  with  that  for  momentum  transfer,  and  the  total  energy  flux  in  the 
turbulent  zone  is 


e- 


PCp« 


ilo 

dy 


For  the  straight  line  velocity  profile 


and 


Therefore 


^ Tg  _  Tpi  —  Toa 
dy  b  . 


(2) 


(3) 


<4) 


(5) 
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and 


q 


PCpf 


(6) 


Prandtl's  expression  for  the  apparent  kinematic  viscosity  is 

6  =  k  (uitiBx  Uii),i2x)b 


(7) 


where  k  is  an  empirical  constant. 


When  Eq.  (7)  is  substituted  into  Eqs.  (5)  and  (6) 

r  =  pk  (uj  -  Ua)  (.8) 

and 

q  -  pcpkdoj  -  To«){uj  -  Ua)  (9) 

As  it  is  desired  to  calculate  shear  stress  and  energy  flux  at  the 
dividing  streamline,  the  flow  properties  must  be  determined  at  that  point 
in  the  mixing  profile.  To  determine  the  dividing  streamline  velocity,  a 
simple  continuity  and  momentum  integral  analysis  was  applied  to  the  two- 
dimensional  mixing  model  of  Fig.  4.  The  results  are  shown  in  Fig.  ^5  for 
equal  stream  temperatures;  the  dividing  streamline  velocity  agrees 
closely  with  the  average  velocity  in  the  mixing  zone.  When  the  effecfof 
stream  temperature  difference  was  treated  for  incompressible  flow,  it 
was  found  that  the  dividing  streamline  velocity  varies  less  than  20  per¬ 
cent  from  the  average  velocity  for  a  range  of  temperature  and  velocity 
ratios  up  to  10.  It  was  therefore  concluded  that,  since  the  separate 
effects  of  compressibility  and  temperature  difference  were  small,  the 
assumption  of  the  dividing  streamline  at  the  midpoint  of  the  mixing  pro¬ 
file  was  justified  for  the  general  case  of  compressible  flow  having  a  stream 
temperature  difference. 


At  the  dividing  streamline  the  velocity,  total  temperature,  and  com¬ 
position  are  averages  of  the  outer  stream  quantities,  and 


Toj  + 


+  c 


Pa 


Uc 


“a 


and  Re  -  -"i  ^ 
^  2 


(The  relationship  for  Rc  is  not  exact  but  is  accurate  for  small  differences 
of  outer  stream  molecular  weights. )  The  density  on  the  dividing  stream¬ 
line,  from  energy  and  state  equations,  is 


Pc 


«PcJ/ 


(10) 
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and  the  shear  ptress  along  the  dividing  streamline  is 

fc  -  Pck  (uj  -  “a) 


(11) 


A  relationship  between  and  qc'  may  be  determined  by  solving 
Eqs.  (5)  and  (6)  for  and  equating 


®Pc  ^ 


(12) 


Thus  Eqs.  (11)  and  (12)  define  the  shear  stress  and  energy  flux  at 
the  dividing  streamline  independent  of  the  mixing  zone  width  and  only  in 
terms  of  the  stream  conditions  outside  the  zone. 


ONE-DIMENSIONAL  FLOW  EQUATIONS 


One- dimensional  equations  for  momentum,  continuity,  and  energy 
may  now  be  written  for  each  stream  over  a  short  interval,  i  x  (Fig.  3), 
and  the  effective  shear  and  heat  flux  are  assumed  to  be  averages  of 
those  quantities  calculated  at  the  interval  end  points. 


When  the  momentum  equation  for  each  stream  is  written  in  terms  of 
the  stream  impulse,  the  sum  of  the  momentum  and  pressure  force  for 
perfect  gases  is 


Faj  =»  Fai  AFj.  —  AFvr  “  Pavg  (Ajj  —  Ajj) 
=  (Ad  “  Ajj)p,  (1  +  yaMaj*) 


where . 


and 


and 


AF„  x=  wDd  AxcJw  Pavg  Ma* 


AF„ 


Fjj  —  Fji  “  AFc  +  Pavg  (Aji  —  Ajj) 


(13) 


(14) 


(15) 


(16) 


=  Aj,  p,  (1  +  yj  Mja’) 
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Since  and  act  on  the  same  surface  areas j  the  total  turbulent 
shear  force  and  total  energy  transfer  rate  acting  over  the  interval  may 
be  related  using  Eq.  (12) 


The  energy  equation  may  then  be  written 

A  Qc  “  Wa  Cp^  (Toa2  "  Toai)  —  Wj  Cpj  (Toji  “  Tojt) 
and  the  weight  flow  equation  for  perfect  gases  becomes 


(17) 


(18) 


and 


(19) 


1  +  yj  ~  ^  Mj/ 
2  ' 


(20) 


By  substitution,  the  preceding  set  of  equations  may  be  reduced  to  six 
unknowns,  p,  i  Aj,,  Man  70821  and  Toja.  with  the  transfer  quantities 
AFc,  AFw»  and  AQo  expressed  in  terms  of  Ax  and  averages  over  the  inter¬ 
val  of  the  six  basic  variables .  The  set  of  equations  could  be  written  for 
solution  with  a  digital  computer;  however,  an  iteration  procedure  amen¬ 
able  to  hand  calculation  has  been  developed. 

The  iteration  procedure  can  be  greatly  facilitated  by  defining  the 
following  parameters.  Thus  Eqs.  (13)  and  (19)  are  combined 


1  +  Xb  M. 


w*  VToi 


M, 


y»  -  1 


‘al 


and  similarly  Eqs.  (16)  and  (19)  are  combined 


(21) 


Fj2  ^  Xj 

Wi  Mj2  yji  + 


(22) 


Equations  (21)  and  (22)  now  present  expressions  for  the  parameter, 
F  /  w  V  To  which  Is  proportional  to  F/F*  of  the  one- dimensional  Fanno 
analysis  for  duct  flow  and  Mach  number.  There  are  two  solutions, 
subsonic  and  supersonic,  for  each  value  of  F  /  w  ,  but  since  the 
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.  r  r 

method  is  concerned  with  a  smooth  transition  of  axial  properties,  the 
Mach  numbers  at  the  end  of  the  calculation  interval  will  be  supersonic 
for  initial  supersonic  conditions  and  subsonic  for  initial  subsonic  con- 
ditipns.  Curves  of  the  supersonic  F  /  w  i/T^  vs  Mach  number  relation¬ 
ship  for  air  and  for  typical  rocket  exhaust  products  are  shown  in  Fig.  6. 
Details  of  the  Iteration  procedure  are  given  In  the  Appendix. 


EMPIRICAL  MIXING  PARAMETERS. 

In  the  development  of  the  simplified  mixing  theory,  only  one  empir¬ 
ical  constant,  k,  was  defined,  other  than  the  wall  skin  friction  coef¬ 
ficient  which  is.  well  known.  An  extensive  literature  survey  revealed 
considerable  discrepancy  in  the  definition  of  mixing  parameters  by 
various  investigators.  Also,  these  empirical  parameters  are  usually 
evaluated  by  correlating  a  particular  theoretical  approach  with  exper¬ 
imental  data,  and  thus  they  reflect  the  peculiarities  of  the  particular 
analysis.  In  addition,  most  of  the  available  data  is  for  the  case  of  con¬ 
stant  pressure  mixing;  however,  it  was  felt  that  turbulent  parameters 
would  not  differ  greatly  for  the  nonconstant  pressure  case  and  that 
information  for  the  constant  pressure  case  would  help  predict  magni¬ 
tudes  and  trends  of  the  empirical  coefficient  for  the  problem  of  inter¬ 
est.  In  any  case,  it  was  necessary  to  determine  proper  values  of  k  for 
the  'simplified  theory  by  correlation  with  experiment. 


Constant  Pressure  Plane  Mixing  with  Zero  Secondary  Velocity 


For  constant  pressure  mixing  between  a  uniform  stream  and  the 
surrounding  ambient  gas,  various  turbulent  exchange  parameters  have 
been  used.  Szablewski  (Ref.  2)  uses  a  similarity  parameter,  a,  which 
by  definition,  is  inversely  proportional  to  the  mixing  zone  spreading 
rate,  dh/d* .  The  various  turbulent  exchange  parameters  are  also  related 
by  (Ref.  2) 

°  ~  V  2k  db/di"  ^23) 

■  Since  is  Inversely  proportional  to  db/dx,  then  k  is  directly  propor¬ 
tional  to  db/dx. 

The  incompressible  value  for  spreading  rate  has  been  established 
by  Tollmien  (Ref.  3)  as  db/dx  =  0.255  for  a  uniform  stream  ihixing  with  a 
surrounding  ambient  gas.  The  value  of  o  12  also  has  been  established 
for  this  case.  Then  Eq.  (23)  yields  k  *  0.0136.  (Szablewski  (Ref.  2)  com¬ 
puted'  the  value  k  «  0.0158 . ) 
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It  is  generally  known  that  increasing  stream  Mach  number  decreases 
the  rate  of  spreading.  Data  obtained  at  Princeton  (Refs.  4  and  5)  and  at 
NACA  (Ref.  6)  are  plotted  in  Fig.  7  for  zero  secondary  flow.  There 
appears  to  be  an  almost  linear  decrease  of  db/dx,  and  therefore  k,  with 
increasing  Mach  number. 

Finito  Secondary  Velocity 


For  finite  secondary  stream  velocity,  ua,  Eq.  (23)  becomes  (Ref.  2) 

=  1 

v  2k  (db/dx)„^=..  Q  (i  -  Ug/ttj  )  (24) 

and  db/dx  varies  as  ^1  -  ua/uj);  k  is  not  affected  by  finite  secondary 
stream  velocities.  Spreading  rate  data  from  Princeton  (Ref.  4)  agree 
well  with  the  predicted  trend  of  db/d*  with  uaAj. 

Effect  of  StToam  Temperature  Difference 

Very  little  experimental  data  have  been  published  concerning  the 
effect  of  stream  temperature  differences  on  turbulent  mixing  rates.  Data 
from  Princeton  (Ref.  4)  for  moderate  temperature  difference  and  sub¬ 
sonic  flow  indicate  that  increasing  the  primary  stream  total  temperature 
tends  to  increase  the  spreading  rate;  however,  insufficient  data  are 
available  to  fully  define  the  trend. 


Constant  Pressure  Mixing  of  on  Axisymmotric  Jot 

For  fully  developed  mixing  of  an  axisymmetric  free  jet  where  there 
is  a  decay  of  centerline  velocity,  somewhat  different  values  of  k  have  been 
observed.  Szablewski  (Ref.  7)  determined  a  value  of  k  «  0.01U5  to  corre¬ 
late  his  theory  with  experimental  data  for  the  diffusion  of  a  fully  developed 
jet.  Warren  (Ref.  8)  correlated  experimental  data  with  his  integral  equa¬ 
tion  theory  for  the  spreading  of  a  fully  developed  axisymmetric  jet  with 
zero  secondary  velocity.  Warren's  incompressible  value  for  k  is  0.0217; 
the  difference  between  this  and  the  value  0.0105  reported  by  Szablewski  re¬ 
flects  differences  in  the  methods  of  analysis.  Warren's  results  for  k  are 
shown  in  Fig.  8  plotted  against  initial  primary  jet  Mach  number  along  with 
the  empirical  relationship  used  to  correlate  the  data. 
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EXPERIMENTAL  INVESTtCATION 

•  ’  i 

EXPERIMENTAL  APPARATUS 
Op*n-Circuit  Tunnel 


The  experimental  program  was  conducted  in  a  small  open- circuit 
wind  tunnel  installed  in  the  Propulsion  Research  Laboratory  <T5C-2)  of 
the  RTF.  A  schematic  of  the  tunnel  is  shown  in  Fig.  1.  Tunnel  air  was 
supplied  from  the  RTF  air  machinery  at  total  pressures  up  to  40  psia. 

A  direct-fired  heater  using  JP-4  was  installed  at  the  inlet  to  the  plenum 
chamber  to  heat  the  tunnel  air  up  to  a  maximum  of  1200°F.  The  flow 
was  accelerated  in  an  8. 1-in.  -diam.  Mach  number  3,  07,  axisymmetric 
nozzle.  A  2.  83 -in.  -diam  centerbody  designed  by  the  method  of  Ref.  9 
to  provide  uniform  annular  flow  past  the  model  base  was  Installed  in  the 
tunnel  nozzle.  The  diffuser  duct  used  for  this  investigation  was  8.  1  in. 
in  diameter  and  had  a  length  of  10  diameters.  The  construction  of  the 
downstream  flange  caused  a  small  (approximately  4  percent)  area  re¬ 
striction  at  the  end  of  the  duct. 

A  small  water-cooled  rocket  motor  utilizing  inhibited  red  fuming 
nitric  acid  (IRFNA)  and  unsymmetrical  dimethyl  hydrazine  (UDMH)  as 
propellants  was  mounted  at  the  downstream  end  of  the  centerbody.  The 
flow  configuration  past  the  end  of  the  centerbody  essentially  simulated 
the  flow  past  a  single  rocket  missile  base  in  free  flight.  A  summary  of 
principal  tunnel  and  rocket  dimensions  and  operating  parameters  is  pre¬ 
sented  in  the  Table. 

Rocket  Instrumentotion 

Rocket  chamber  pressure  and  injector  head  pressures  were  meas¬ 
ured  by  strain-gage  transducers  and  recorded  on  a  recording  oscillo¬ 
graph.  Fuel  and  oxidizer  flow  rates  were  measured  by  turbine-type 
flowmeters  and  recorded  on  the  oscillograph.  The  sequence  of  propel¬ 
lant  valve  operation  was  also  recorded  on  the  oscillograph. 

Rocket  chamber  pressure  and  propellant  system  tank  pressures  were 
measured  by  strain- gage  transducers  and  recorded  on  strip  charts  for 
inspection  during  rocket  operation.  Fuel,  oxidizer,  and  cooling  water 
flow  rates  were  also  monitored  visually  by  digital  time  and  frequency 
meters. 


Aerodynamic  Instrumentation 

Pressures  necessary  for  setting  tunnel  operating  conditions  (tunnel 
total  pressure,  RTF  exhauster  pressure,  spray  section  pressure,  and 
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nozzle  exit  static  pressure)  were  measured  on  high  accuracy  diaphragm 
gages.  These  pressures,  as  well  as  the  static  pressure  distribution 
from  the  nozzle  exit  to  the  spray  section,  were  also  measured  by  mer¬ 
cury  manometers  referenced  to  atmosphere  and  recorded  photographi¬ 
cally.  The  estimated  error  of  the  pressures  measured  on  the  manometers 
is  ±0,  05  psi. 

Mach  number  distribution  over  the  nozzle  annulus  was  determined 
during  initial  calibration  runs  by  a  5 -station  total  pressure  rake  installed 
just  upstream  of  the  end  of  the  centerbody.  The  pressures  were  meas¬ 
ured  with  mercury  manometers  and  recorded  photographically. 

Tunnel  total  temperature  was  measured  by  iron- cons tantan  thermo¬ 
couples  read  out  visually  on  temperature -calibrated  galvanometers. 


EXPERIMENTAL  PROCEDURE 

The  following  procedure  was  used  to  obtain  the  data  on  supersonic 
two- stream  mixing: 

1.  Tunnel  flow  with  a  very  low  back  pressure  was  established  so 
that  supersonic  flow  existed  in  the  entire  straight  section  of  the 
diffuser. 

2.  The  rocket  fuel  and  oxidizer  valves  were  opened,  and  ignition 
was  accomplished  when  the  hypergolic  propellants  were  injected 
into  the  chamber. 

3.  The  rocket  firing  continued  for  20  to  30  sec  to  allow  stabilization 
of  the  instrumentation.  Pictures  of  the  manometers  were  taken 
automatically  every  three  seconds  during  the  rocket  firing  to 
obtain  a  record  of  the  manometer  response. 


DISCUSSION  OF  RESULTS 


To  verify  that  the  simplified  theory  gives  an  adequate  explanation  of 
the  bounded  two-stream  mixing  process,  correlations  were  made  with 
Mikhail's  experimental  data  (Ref.  1)  for  incompressible  flow.  (A  sche¬ 
matic  of  Mikhail's  experimental  apparatus  is  shown  in  Fig.  9.)  The  simpli¬ 
fied  theory  was  also  correlated  with  Mikhail's  more  elaborate  integral 
equation  theory  for  incompressible  mixing.  The  one-dimensional  flow 
equations  were  written  in  a  form  applicable  to  incompressible  flow,  and 
calculations  were  carried  out  using  various  values  of  the  turbulent  ex¬ 
change  parameter,  k,  for  the  initial  conditions  shown  in  Fig.  9.  It  was 
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found  that  the  best  correlation  was  obtained  using  k  =  0.010,  which  is 
essentially  the  value  used  by  Szablewski  (Ref.  7)  for  the  constant  pres¬ 
sure  diffusion  of  a  fuHy  developed  incompressible  jet. 

A  comparison  of  Mikhail's  experimental  axial  pressure  distribution 
with  theoretical  distributions  calculated  from  the  simplified  theory  and 
from  Mikhail’s  theory  is  shown  in  Fig,  10.  Both  theoretical  pressure 
distribution  curves  differ  considerable  from  the  experimental  curve  for 
the  first  eight  duct  diameters.  Mikhail  explains  this  discrepancy  as 
caused  by  a  low  pressure  region  in  the  center  of  the  primary  jet  which 
extends  about  the  first  eight  duct  diameters ;  however,  this  explanation 
violates  the  assumption  in  both  theories  of  constant  pressure  across  any 
section.  Other  investigators  have  also  found  the  same  shape  for  the 
experimental  pressure  distribution  curve,  and  it  appears  typical  for 
low  velocity  jet  pump  configurations. 

A  comparison  of  the  experimental  centerline  velocity  distribution 
with  theoretical  distributions  calculated  using  the  two  theories  is  shown 
in  Fig.  11.  The  simplified  theory  correlates  the  experimental  data 
poorly  for  the  first  six  to  eight  duct  diameters  but  shows  acceptable 
agreement  farther  downstream.  Mikhail's  theory  gives  excellent  corre¬ 
lation  with  experimental  velocity  distribution  data  since  each  of  his 
three  empirical  mixing  parameters  was  determined  by  matching  these 
theoretical  and  experimental  centerline  velocity  distributions.  Although 
the  simplified  theory  presented  in  this  report  only  approximately  pre¬ 
dicts  experimental  velocity  and  pressure  distributions  for  the  first  six 
to  eight  duct  diameters,  the  behavior  farther  downstream  is  predicted 
closely  and  agrees  well  with  the  predictions  of  Mikhail's  more  elaborate 
theory.  It  can  be  concluded  that  the  simplified  theory  gives  a  resonably 
accurate  overall  picture  of  the  mixing  process  and  may  be  used  to  opti¬ 
mize  low  velocity  mixing  tube  configurations. 


COMPRESSIBLE  TWO-STREAM  FLOW 

A  basic  assumption  of  the  simplified  analysis  is  that  the  two  streams 
are  initially  uniform  and  parallel.  However,  the  actual  wind  tunnel  con¬ 
figuration  used  was  not  designed  specifically  to  study  the  mixing  problem 
and  violates  this  assumption  to  a  certain  extent.  A  schematic  of  the 
model  base  geometry  and  flow  just  downstream  of  the  rocket  is  shown  in 
Fig.  12.  To  apply  the  simplified  mixing  theory  to  the  initial  flow  con¬ 
ditions  of  the  actual  tunnel,  it  was  assumed  that  at  some  section  down¬ 
stream  the  two  streams  become  uniform  and  parallel.  The  flow  condi¬ 
tions  at  this  hypothetical  section  were  calculated  assuming  that  each 
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Stream  undergoes  an  isentropic  area  change  to  equalize  pressures  in  the 
two  streams.  These  flow  conditions  are  than  used  as  the  initial  condi¬ 
tions  for  the  mixing  analysis. 

Axial  duct  static  pressure  distributions  for  the  Mach  number  3..  07 
tunnel  operating  at  a  tunnel  total  temperature  of  1100*R  are  shown  in 
Figs.  13a,  b,  and  c.  The  experimental  rocket- off  distributions  shown 
indicate  that  the  tunnel  back  pressure  was  sufficiently  low  to  allow  super¬ 
sonic  flow  over  the  entire  duct.  In  the  case  of  rocket- on  operation,  the 
high  initial  duct  static  pressure  just  downstream  of  the  model  was 
caused  by  the  impingement  and  reflection  at  the  wall  of  the  shock  formed 
at  the  intersection  of  the  two  streams.  After  about  two  to  three  duct 
diameters,  the  magnitude  of  the  disturbances  was  considerably  de¬ 
creased,  and  a  reasonable  average  static  pressure  distribution  may  be 
drawn  through  the  experimental  points.  The  initial  shock  increased  in 
strength  with  lower  tunnel  static  pressure  (Increased  underexpansion  of 
the  rocket  nozzle),  and  the  assumption  of  uniform  initial  stream  con¬ 
ditions  was  violated  to  a  greater  extent  as  tunnel  pressure  decreased. 

When  the  simplified  mixing  theory  was  applied  to  the  experimental 
data  (Fig.  13),  it  was  found  that  the  best  correlation  was  obtained  for 
k  ^  0.010,  the  value  used  for  incompressible  flow.  Apparently  the  trend 
of  decreasing  k  with  increasing  Mach  number  is  equally  opposed  by  the 
trend  of  increasing  fc  with  increasing  stream  temperature  difference  for 
the  particular  experimental  parameters. 

Experiments  were  also  made  with  a  lower  tunnel  total  temperature 
(500-600"R).  In  this  case  thermal  choking  occurred  in  the  10-diam  duct, 
even  though  the  back  pressure  was  maintained  at  a  value  lower  than  the 
duct  exit  pressure.  The  choking  caused  the  tunnel  nozzle  flow  to  break 
down,  and  the  resulting  static  pressure  distributions  were  indistin¬ 
guishable  from  those  which  occurred  when  the  mixing  of  the  streams  was 
induced  by  a  fully  developed  duct  shock  system.  The  simplified  theory 
was  applied  to  the  cases  where  choking  was  observed,  and  the  resulting 
outer  stream  axial  Mach  number  distributions  are  shown  in  Fig.  14,  At 
the  10-diam  station,  theoretical  values  for  Ma  are  1.25  and  1.48,  re¬ 
spectively,  for  tunnel  pressures  of  19.  9  and  37.  8  psia.  The  4-percent 
area  reduction  at  the  downstream  flange,  however,  would  cause  choking 
to  occur  at  approach  Mach  numbers  of  1.2  or  1.3  rather  than  1.  0.  Also, 
the  theoretical  Mach  number,  Ma,  is  based  on  the  discontinuous  velocity 
profile  in  Fig.  3  and  represents  an  average  Mach  number  for  the  outer 
airstream.  The  total  pressure  loss  caused  by  the  shocks  formed  at  the 
intersection  of  the  initially  non-uniform  streams  would  also  cause  choking 
to  occur  at  theoretical  values  of  Ma  somewhat  greater  than  1.0.  Theoret¬ 
ical  Mach  number  distributions  for  1100"R  tunnel  total  temperature  where 
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no  choking  occurred  are  shown  in  Fig.  15;  in  this  case  the  theoretical 
value  of  Ma  is  approximately  2.  0  at  the  lO'diam  station.  It  can  thus  be 
concluded  that  the  simplified  theory  indicates  the  correct  trends  of  outer 
stream  Mach  number  and  axial  static  pressure  and  that  the  essential 
features  of  the  mixing  process  are  preserved  in  the  analysis. 


CONCLUSIONS 


The  following  conclusions  may  be  drawn  from  the  portion  of  the 
two- stream  diffuser  investigation  reported  herein; 

1.  The  simplified  treatment  of  constant  area  two-stream  mixing 
preserves  the  essential  features  of  the  process  and  correlates 
well  with  other  approaches  for  the  case  of  low  velocity  flow. 

2.  The  simplified  theory  explains  the  general  behavior,  with  re¬ 
spect  to  Mach  number  and  static  pressure,  of  the  constant  area 
mixing  of  two  supersonic  streams. 

3.  Thermal  choking  in  a  constant  area  duct  can  be  caused  by  tur¬ 
bulent  mixing  between  a  supersonic  rocket  stream  and  a  sur¬ 
rounding  supersonic  airstream,  even  when  the  duct  back 
pressure  is  maintained  at  a  low  value. 
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APPENDIX 

ITERATION  PROCEDURE  FOR  MIXING  ANALYSIS 


For  the  general  case  of  bounded  mixing  of  a  rocket  stream  with  a 
surrounding  airstream,  the  following  iteration  procedure  was  used: 

A.  First  Iteration 

All  parameters  at  the  beginning  of  the  calculation  interval  are  given. 
Choose  AQc  which  defines  Toa»  and  Toj*  from  Eq.  (18).  As  a  first  approx¬ 
imation  assume  that  the  average  transfer  quantities  over  the  interval  are 
those  calculated  at  the  initial  section. 


1. 


2, 

3. 


4. 


5. 


Calculate  rei  from  Eq,  (11)  using  proper  value  of  the  empirical 
parameter  k. 


Use 


~  um 

Tgjl  ~  Tggl 


for 


T  -  “  T 

oj  ®^avg 


and  calculate  A  Pc  from  Eq.  (17). 


Use  rci  for  Dji  for  and  solve  for  Ax  from  Eq.  (15). 

Calculate  AF*  from  Eq.  (14)  using  initial  stream  conditions  and 
Ax  calculated  in  step  3. 

Neglect  the  pAA  force  and  calculate  Fa,  and  Fj,  from  Eqs.  (13) 
and  (16). 


6. 


7. 


8. 


Evaluate  the  parameters  Fgj  /  Wg  V'l'oa,  and  Fj,  /  wj  VToji 

and  determine  Ma,  and  Mj,  from  a  plot  similar  to  Fig.  6  for  the 

appropriate  gas  properties. 


Calculate  na,  and  uj, ,  and  determine 


for  the  interval. 


Solve  Eqs.  (13)  and  (16)  simultaneously  for  p,  and  Aj, . 


9.  Calculate  pc,  and  evaluate  ^c,  from  Eq.  (11) 


B.  Second  Iteration 

1.  Recalculate  AFc  from  Eq.  (17)  using  average  quantities  found  in 
the  first  iteration. 


2.  Solve  for  Ax  from  Eq.  (15),  again  using  average  quantities  over 
the  interval. 


3.  Calculate  AFw  from  Eq.  (14)  using  average  outer  stream  condi¬ 
tions  and  Ax  calculated  in  step  2. 

4.  Calculate  Pavg  force  from  the  first  iteration  and  evaluate  Fa, 
and  Fja . 
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5.  Repeat  steps  (6)  to  (8)  of  first  iteration. 

6.  Solve  for  Ax  from  Eq.  (15). 

The  iteration  procedure  is  repeated  until  the  values  for  and  Fj, 
calculated  at  the  end  of  an  iteration  agree>  to  the  desired  accuracy,  with 
those  used  at  the  beginning  of  that  iteration. 

After  the  flow  properties  at  the  end  of  the  initial  interval  are  cal¬ 
culated  to  the  desired  accuracy,  the  process  is  repeated,  using  the  end 
conditions  of  the  first  interval  as  initial  conditions  for  the  next.  It  was 
found,  for  the  parameters  of  the  experimental  apparatus,  that  incre¬ 
ments  of  Toa  of  lOO^R  gave  adequate  accuracy. 

If  any  simplifying,  condition  exists,  such  as  no  wall  friction,  con¬ 
stant  stream  temperature,  or  incompressible  flow,  the  iteration  pro¬ 
cedure  is  greatly  simplified.  It  is  also  possible  to  include  duct  area 
change;  however,  the  iteration  procedure  is  more  complicated. 
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TABLE 

TUNNEL  DIMENSIONS  AND  OPERATING  PARAMETERS 


Tunnel  Nozzle  Exit  Diameter 

Centerbody  Diameter  (nozzle  exit  plane) 

Rocket  Nozzle  Exit  Diameter 

Rocket  Nozzle  Area  Ratio 

Rocket  Nozzle  Half  Angle  (conical  nozzle) 

Nominal  Rocket  Chamber  Pressure 

Nominal  Rocket  Weight  Flow 

Nominal  Rocket  Oxidizer- Fuel  Ratio 

Nominal  Rocket  Vacuum  Thrust  (calculated) 

Theoretical  Rocket  Chamber  Temperature 

Theoretical  Rocket  Exhaust  Gas  Specific  Heat  ( cpj ) 

Theoretical  Rocket  Exhaust  Gas  Molecular  Weight 

Theoretical  Rocket  Exhaust  Gas  Ratio  of  Specific 
Heats  (y) 


8. 1  in. 

2.  83  in. 

1.  25  in. 

8 

15  deg 
520  psia 

0.  52  Ibm/sec 
2.8 
149  Ibf 
5500®R 

0.  46  Btu/lbjn"R 
23.  3 
1.  23 
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TUNNEL  AIR  FROM 
DIRECT-FIRED  PREHEATER 
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Fig.  1  Schematic  of  Tunnel  for  Two-Stream  Diffuser  Investigation 
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Fig*  4  Mixing  Model  for  Calculation  of  Local  Rates  of  Heat  ond  Momentum  Transfer 
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Figt  5  Dividing  Streamline  Velocity  Ratio 


w/T( 


AEDC-TR-61-18 


MACH  NUMBER 

Fig.  6  Mach  Number  Funetien  far  Air  and  Typical  Rocket  Preducte 
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MACH  NUMBER 

Fig.  7  Variation  of  Mixing  Zono  Spreading  Rate  with  Mach  Number  for  Zero  Secondary  Velocity 


33 


TURBULENT  MIXING  PARAMETER 


AEDC-TR-6».18 


34 


9  Schematic  of  MikhaiTs  Experimentol  Apparatus 
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DIMENSIONLESS  PRESSURE  RISE, 
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Fig.  10  Comparison  of  Theory  wilh  Mikhail's  Experimental  Dote 
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Fig.  11  Comparison  of  Thaoreticol  and  Experrmentol  Centerline  Velocity  Distribution 
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Fig.  7  Continued 


DUCT  STATIC  PRESSURE 
NOZZLE  EXIT  STATIC  PRESSURE 


Fig.  7  Concluded 
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Fig.  14  Theoretical  Presiure  and  Mach  Number  Distributions  for  Mach  Number  3.07  Tuntwl 
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Fig-'  15  Thaefctleol  Outer  Stream  Mach  Number  Distribution  for  Moeh  Number  3,07  Tunnel 
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